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(Dated: August 29, 2008) 

We report the observation of the doubly strange b baryon in the decay channel Q'^ ^ J/^ Q~ , 
with J/ip —> fi^ n~ and fl~ AK~ —> {p-K~)K~ , in pp collisions at y/s = 1.96 TeV. Using 
approximately 1.3 fb~^ of data collected with the DO detector at the Fermilab Tevatron Collider, we 
observe 17.8±4.9 (stat.)±0.8 (syst.) 0.^ signal events at a mass of 6.165±0.010 (stat.)±0.013 (syst.) 
GeV. The significance of the observed signal is 5.4fT, corresponding to a probability of 6.7 x 10~* of 
it arising from a background fluctuation. 
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PACS numbers: 14.20.-c, 14.20.Mr, 14.65.Fy 



The n~ baryon, composed of three strange quarks, 
played an important historical role in our understanding 
of the basic structure of matter. Its discovery in 1964 [1| 
at a mass predicted from SU(3) symmetry breaking was 
a great success for the theory The fl^ (bss) (charge 
conjugate states are assumed throughout this Letter) is 
a predicted heavy cousin of the rt~ with a b quark re- 
placing one of the three strange quarks. While the fi" 
has = 3/2+, the ground state fl^ is expected to have 
= 1/2+, a mass between 5.94 — 6.12 GeV and a life- 
time such that 0.55 < T{n-)/T{B°) < 1.10 

In this Letter, we report the first observation of the f2(^ 
baryon, fully reconstructed from its decay fl^ — > J/ijj , 
with J/^p — > /i+M^, il" AK~ and A pn^ . The 
analysis is based on a data sample of 1.3 fb^^ collected 
in pp^ collisions at -y/s = 1.96 TeV with the DO detec- 
tor [J] at the Fermilab Tevatron Collider. The detector 
components most relevant to this analysis are the cen- 
tral tracking system and the muon spectrometer. The 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT) inside a 
2 Tesla superconducting solenoid. The SMT is optimized 
for tracking and vertexing over the pseudorapidity region 
I77I < 3 while the CFT has coverage for \ri\ < 2. A liquid 
argon and uranium calorimeter provides coverage up to 
I77I < 4.2. The muon spectrometer covers \ri\ < 2. 

The J/^Q- J/^KK- J/ijpn- K- de- 

cay topology is similar to that of the J/V'S^ 
J/ip Att~ —^ J/ippn^ TT^ decay with in place of 
and in place of tt^. Consequently, the reconstruc- 
tion of the J/ip and A and their selection discussed below 
follow closely the analysis that led to the first direct ob- 
servation of the baryon However, in this analysis 
we use a multivariate selection for the owing to the 
smaller signal to background ratio compared to that for 
the S~ in the analysis. We use the pythia Monte 
Carlo (MC) program 6] to generate fi^ and the evtgen 
program [71] to simulate il^ decays. The mass and 
lifetime are set to be 6.052 GeV and 1.54 ps respectively. 
The generated events are subjected to a GEANT [1| based 
DO detector simulation, and to the same reconstruction 
and selection programs as the data. We optimize the fi" 
selection using MC Sl^ events for the signal and a sam- 
ple of J/ip{AK'^) data (referred to below as wrong-sign 
events) for the background, while keeping the J/ijjQ,^ 
data blinded. Once all selection criteria have been deter- 
mined, we apply them to the J /ip data. 

We begin the event selection by reconstructing J /ip 
/i+/i~ candidates from two oppositely charged muons 
with transverse momentum (pr) greater than 1.5 GeV 
that are compatible with being from a common vertex. 
Muons are identified by tracks reconstructed in the cen- 
tral tracking system that are matched with either track 



segments in the muon spectrometer or calorimeter en- 
ergy deposits consistent with a minimum ionizing parti- 
cle. Events must have a well-reconstructed pp interac- 
tion point that we take to be the JljJ" production ver- 
tex and a J/V' IJ-^ candidate in the mass window 
2.75 < < 3.40 GeV. Events with J/^p candidates are 
reprocessed with a version of the track reconstruction al- 
gorithm that increases the efficiency for tracks with low 
Pt and high impact parameters. 

We form A — > p'K~ candidates from two oppositely 
charged particles, each withpr > 0.2 GeV, that are con- 
sistent with having originated from a common vertex. 
The two tracks are required to have a total of no more 
than two hits in the tracking system before the recon- 
structed p'K~ vertex. The impact parameter significance 
(the impact parameter with respect to the pp interaction 
point divided by its uncertainty) must exceed four for at 
least one of the tracks and three for the other. The track 
with the higher pT is assumed to be the proton. MC 
studies show that this assignment leads to the correct 
combination nearly 100% of the time. Furthermore, we 
require the A transverse decay length to be greater than 
four times its uncertainty and the proper decay length 
to exceed ten times its uncertainty, where the transverse 
decay length is the distance between the production and 
decay vertices in the transverse plane while the proper de- 
cay length is the transverse decay length corrected by the 
Lorentz boost calculated from pt(A). A candidates must 
have reconstructed masses between 1.108 and 1.126 GeV. 




M(AK) (GeV) M(AK) (GeV) 



FIG. 1: The invariant mass distribution of the KK pair before 
(a) and after (b) the BDT selection. Filled circles are from the 
right-sign KK~ events while the histogram is from the wrong- 
sign KK'^ events without any additional normalization. 

We combine A candidates with negatively charged par- 
ticles (assumed to be kaons) to form fi~ AK~ decay 
candidates. The A and the kaon are required to have 
a common vertex. The $7" candidates must have trans- 
verse decay length significances greater than four and the 
uncertainties of the proper decay lengths less than 0.5 cm. 
These two requirements reduce backgrounds from com- 
binatorics and mismeasured tracks. The S~ baryon has 
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a mass of 1.322 GeV [9| and decays into Att~. If the 
kaon mass is assigned to the pion, this decay could be 
a major background for AK^ . To ehminate this 

background, we remove candidates with Att" mass less 
than 1.34 GeV. Figure [Ija) shows the mass distribution 
of the reconstructed fl' AK^ candidates after these 
selections. The distribution of wrong-sign AK~^ events is 
also shown. An excess of events above the background 
around the expected fl~ mass of 1.672 GeV Q is visible 
in the distribution of the right-sign AK~ events. 

To further enhance the Q~ signal over the combina- 
torial background, kinematic variables associated with 
daughter particle momenta, vertices, and track qualities 
are combined using Boosted Decision Trees (BDT) |l3, 
ll| . The AK~ mass distribution after the BDT selection 
is shown in Fig. W[h). The BDT selection retains 87% 
of the ri~ signal while rejecting 89% of the background. 
The enhanced 17 ~ mass peak is evident in the distribu- 
tion. A AK~ pair is considered to be a n~ candidate if 
its mass is in the range of 1.662 — 1.682 GeV. 

To select J/tp il^ candidates, we develop selec- 

tion criteria using the MC events as the signal and 
the data wrong-sign events as the background. The back- 
ground events are formed by combining J/ip candidates 
with AK+ pairs with mass between 1.662 and 1.682 GeV. 
We form fi^ — > J/tpn~ decay candidates from J/ip and 
Cl^ pairs that are consistent with being from a common 
vertex. We require the uncertainty of the fi^ proper de- 
cay length to be less than 0.03 cm and impose a minimum 
Pt cut of 6 GeV on the fl^ candidates. Finally, J/ip and 
ri~ daughters from the decays are expected to be 
boosted in the direction of the 17^; therefore, we require 
the opening angle in the transverse plane between the 
J/ip and the fi" to be less than 7r/2. 

We then apply the above selections to the right-sign 
events in the data to search for the baryon in the 
mass window between 5.6 and 7.0 GeV. This range is 
chosen since 5.624 GeV is the mass of the lightest b 
baryon, the Ah, and the upper limit of 7.0 GeV is nearly 
1 GeV higher than the predicted mass [5] . We calcu- 
late the candidate mass using the formula M(f7^) — 
M{J/i}fl')-M{^i+fi-)-M{AK-) + M{J/i}) + M{n-). 
Here MiJ/ipn'), M{fj,+^i,-), and M{AK-) are the re- 
constructed masses while M{J/tp) and M{n~) are taken 
from Ref. Q . This calculation improves the mass resolu- 
tion of the MC events from 0.080 GeV to 0.034 GeV. 
In the mass search window, we observe 79 candidates in 
the data with the mass distribution shown in Fig. [Ifa). 
An excess of events near 6.2 GeV is apparent. No such 
structure, however, is seen in the corresponding mass dis- 
tribution (Fig.[2Ub)) of the 30 wrong-sign events. 

Assuming the excess is due to the production, we fit 
fif^ candidate masses with the hypothesis of a Gaussian 
signal plus a flat background using an unbinned likeli- 
hood method. We fix the Gaussian width to 0.034 GeV, 
the width of the MC fl^ signal. The fit gives an il^ mass 



of 6.165±0.010 (stat.) GeV and a yield of 17.8±4.9 (stat.) 
signal events. To assess the significance of the excess, 
we first determine the likelihood Cs+t of the signal plus 
background fit above and then repeat the fit with only 
the background contribution to find a new likelihood Cb- 
The logarithmic likelihood ratio ^2 \n{£s+b/^b) yields a 
statistical significance of 5.4cr, equivalent to a probability 
of 6.7 X 10~* that the background could fluctuate with 
a significance equal to or greater than what is observed. 
Fitted yields for positively and negatively charged can- 
didates are 6.2 ±3.1 (stat.) and 12.0 ± 3.9 (stat.) n~ , 
respectively. 
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FIG. 2: 



distribution of the candidates 



(a) The M(n^ 

after all selection criteria. The dotted curve is an unbinned 
likelihood fit to the model of a constant background plus a 
Gaussian signal. The mass distributions for the wrong-sign 
background (b), the Jl" sideband events (c), and the A side- 
band events (d). 

Various checks have been performed to ensure that 
the observed resonance is genuine. (1) We apply 
the event selection to data events in the sidebands 
of the reconstructed fi" and A resonances separately. 
The J/i/j (p7r~) K~ mass distributions of these sideband 
events are shown in Figs, [^c) and[2)^d). No significant 
excess is present in either distribution. (2) We inves- 
tigate the possibility of a false signal due to residual b 
hadron backgrounds by applying the final Q.^^ selection 
to MC B- J/ijjK*- J/iPK^TT-, J/V'S", 
and Ab ~* J/ip A samples with equivalent integrated lu- 
minosities significantly greater than that of the analyzed 
data. No indication of any resonance is observed in the 
reconstructed J/^^l" mass distribution. (3) We check 
the mass distributions of the fi^ decay products. For 
il^ candidates within a ±3(7 mass window around the 
observed peak, we relax the mass requirements on the 
H.^ and A candidates and perform a fit to each mass dis- 
tribution. The numbers of the f2~ and A candidates from 
the fits are consistent with the observed number of fl^ 
signal events. (4) We replace the BDT selection with in- 
dividual cuts on the most important variables according 
to the BDT optimization and confirm the existence of a 
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peak with a comparable event yield but a higher back- 
ground at a mass consistent with that observed using the 
BDT. (5) We test the robustness of the peak by varying 
selections such as the 5~ veto, A and n~ mass windows, 
A transverse decay requirements, BDT selection, and the 
requirement on pt{^^)- All the above studies confirm 
the existence of the resonance. 

Potential sources of systematic uncertainties on the 
measured fi^ mass include event selection, signal and 
background models, and momentum scale. Varying the 
selection criteria and applying a set of cuts on individ- 
ual kinematic variables lead to a maximum change of 
0.012 GeV in the fitted mass. Using a linear function as 
the background model results in negligible change in the 
mass. Varying the Gaussian width in the signal model 
between 0.028 and 0.040 GeV changes the fitted mass 
by at most 0.003 GeV. When a tighter selection is ap- 
plied to enhance signal over background, we can float 
the width of the signal model in the fit. This leads to 
a mass shift of 0.002 GeV and a fitted signal width of 
0.033 ± 0.010 GeV, consistent with the MC expecta- 
tion. To study the effect of the track momentum scale 
uncertainty on the measured fi^ mass, we reconstruct 
the higher statistics Ab — > J/tp A decays and measure 
the Ab mass for different minimum px requirements on 
the Ab daughter particles. We compare these measure- 
ments to the world average value of the Ab mass [9] and 
take the maximum deviation of 0.004 GeV as a system- 
atic uncertainty. Adding in quadrature, we get a total 
systematic uncertainty of 0.013 GeV to obtain a mea- 
sured r^^mass: 6.165 ± 0.010 (stat.) ± 0.013 (syst.) GeV. 
Similarly, we estimate the systematic uncertainty on the 
flf^ yield by varying the signal and background mod- 
els in the fit. The observed maximum change of 0.8 
is assigned as the systematic uncertainty on the yield: 
17.8 ± 4.9 (stat.) ± 0.8 (syst.). In aU these studies, the 
signal significance remains above 5a. 



Data 

MC signal + data bkgd 




0.1 0.2 0.3 0.4 

Proper decay length (cm) 

FIG. 3: The distribution of the proper decay length of the 
candidates in the ±3(t mass window around the observed 
peak along with the expected distribution from the MC 
signal with a lifetime 1.54 ps plus the data background events. 



shown. The background distribution is modeled using 
events in the fi^ sidebands of 5.8 — 6.0 and 6.4 — 6.6 GeV. 
Despite the low statistics, the data distribution contains 
significantly more positive than negative decay lengths as 
expected and consistent with a weakly decaying b hadron. 

We calculate the fi^ production rate relative to that 
of the 0. The selection efficiency ratio e(ri^ 
J/^n-)/e{E^ J/ipE-) is found to be 1.5±0.2 (stat.) 
assuming inclusive and decays. The higher effi- 
ciency for the is due primarily to a harder px spec- 
trum of the kaon from the Q~ decay than that of the 
pion from the decay and a shorter lifetime of the f2~ 
compared to the S~. By using the reported events [3| 
and the observed fi^ yield here, we estimate 



n 



/(fc->»b') Brjn- ^ j/^jn-) 



to be 7^ = 0.80±0.32 (stat.)I^:22 (syst.). Here f{b n^) 
and f{b S^) are the fractions of b quarks that 
hadronize to form and S^, respectively. The sys- 
tematic uncertainty includes contributions from the sig- 
nal yields as well as the efficiency ratio. Using r(r2^ — + 



J/ipn-)/T{Ef; J/ipE-) = 9.8 []J|, the central val- 
ues of t(S^) = 1.421;q;24 ps 0, the TZ value above, 
and T(ri^) in the range of 0.83 — 1.67 ps 3, we obtain 
f{b ^ n-)/f{b ^ E-) « 0.07 - 0.14. 

In summary, by analyzing 1.3 fb ^ of data collected by 
the DO experiment in pp collisions at ^/s — 1.96 TeV at 
the Fermilab Tevatron Collider, we have made the first 
observation of the doubly strange b baryon in the fully 
reconstructed decay mode — s- J/'il)Q~ with J/iJj — > 
/i"*"/!^, fl^ AK^ and A pn^ . We measure the fi^ 
mass to be 6.165 ± 0.010 (stat.) ± 0.013 (syst.) GeV. The 
significance of the observed signal is greater than 5cr. 
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Figure [3] shows the distribution of the proper decay 
length of the il'^ candidates observed in the ±3tT mass 
window around the fitted fi^ mass. The distribution of 
the MC rijj' signal plus the data background events is also 
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